SUBSTITUTE SPECIFICATION (Marked-Up) 

ArMETHODS OF ESTIMATING PRECIPITATION CHARACTERISTICS 

Related Application 

[0001] This is a §371 of International Application No. PCT/FR2004/002692, with an inter- 
national filing date of October 21, 2004 (WO 2005/040853 Al, published May 6, 2005), which 
is based on French Patent Application No. 03/12299, filed October 21, 2003. 

Technical Field 

[0002] Th e pr e s e nt This invention relates to ar-methods of estimating precipitation char- 
acteristics and a in one p articula r aspect, precipitation rate for solid precipitation. 



Background 

[0003] It is known from th e stat e of th e art that the characteristics of a radar image can be 
used to estimate rainfall characteristics. In particular, Europ e an Pat e nt EP 1 049 944 describes a 
technique for estimating rainfall using a radar. In that patent, provision is made for the following 
steps: 

[[-]] using thata dual-polarized radar to measure the differential phase shift (Odp) and 
the apparent reflectivity Z in at least one of the H (horizontal) or V (vertical) polariza- 
tions and over a given range {{ro, ri}) of path radii relative to saklthe radar; 
[[-]] determining an estimate of the value No* representative of the size distribution of 
the rain drops, on the basis of the differential phase shift difference over the range ro to 



ri, and on the basis of an integral of a function of the apparent reflectivity Z, over the 
range f{ro, ri}); and 

[[-]] deducing the value of the precipitation rate at a point on the basis of No* and of 

the apparent reflectivity at sakfche point. 
[0004] A method of estimating rainfall is also known from PCT Pat e nt WO_03/007016, 
which describes a method of estimating a precipitation rate by means of a dual-polarized radar, 
that method being characterized by the following steps: 

[[-]] using thata dual-polarized radar to measure the differential phase shift Odp and 

the-attenuated reflectivity Z in at least one of the H or V polarizations and over a given 

range £(rl, r0}) of path radii r relative to smdthe radar; 

[[-]] determining an estimate of the value K(r0) of the specific attenuation at rO on the 
basis of the profile of the attenuated reflectivity measured in this way, and on the basis of 
the differential phase shift difference over the range rO to rl; determining an estimate 
K(r) of the specific attenuation at r as a function of the attenuation K(r0) determined in 
this way and of the attenuated reflectivity profile Z(r); and determining the precipitation 
rate R(r) once K(r) is known. 
[0005] Those various solutions make it possible to characterize liquid precipitation and te 
estimate precisely rainfall rate (in millimeters per hour (mm/h), but they do not make it possible 
to estimate the characteristics of solid precipitation such as snow. 

Summary 

[0006] This invention relates to a method of estimating precipitation characteristics including 
acquiring a radar image including at least a vertical plane of a precipitation zone, processing a 
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vertical profile to generate digital signals representative of reflectivity in vertical direction z. 
integrating the signals representative of reflectivity by assimilation of a reflectivity vertical 
profile in an aggregation model to generate a signal representative of the profile in the vertical 
plane of a mean particle diameter weighted by mass of each particle, and determining concen- 
tration of the solid particles on the basis of signals previously determined. 

[0007] An obj e ct of th e inv e ntion is to provid e a nov e l m e thod that r e m e dies that drawback. 

[0008] To this e nd, and in its most g e neral acc e ptation, th e inv e ntion provid e s a m e thod of 

e stimating pr e cipitation charact e ristics and in particular pr e cipitation rat e for solid precipitation, 
said m e thod comprising an acquisition and proc e ssing step consisting in acquiring a radar imag e 
including at l e ast a v e rtical plan e of a pr e cipitation zon e , and in proc e ssing a v e rtical profil e so as 
to d e liv e r digital signals r e pr e s e ntativ e of refl e ctivity in th e v e rtical dir e ction h, said m e thod 
being charact e rized in that it further compris e s an int e gration st e p consisting in int e grating said 
signals r e pr e s e ntativ e of r e fl e ctivity so as to d e liv e r a signal repr e s e ntativ e of th e profil e in the 
vertical plan e of a th e m e an particl e diam e t e r w e ight e d by th e mass of e ach particle, and a det e r 
mination st e p consisting in d e termining th e conc e ntration of th e solid particl e s on the basis of th e 
signals comput e d in th e pr e c e ding st e ps. 

{0009] Pr e f e rably, th e int e gration st e p consists in d e t e rmining th e variabl e Z(h) of th e radar 

observabl e in mm ^ j as a function of the altitud e h on the basis of said radar imag e , and in 
d e t e rmining said m e an diameter Dm(h) of th e particl e s by resolving th e following equation: 
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wh e r e : 



Z is th e radar obs e rvable to be inv e rt e d in mm m f 



S ffl is in m e t e rs (m); 
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a and b ar e co e ffici e nts specific to particles of tho "aggregat e " type; for example, tho co 
e ffici e nt a is equal to 3518 4 and th e co e fficient b is equal to 3.16; 

&<ff is th e co e ffici e nt of e ff e ctiv e ness of tho aggregation proc e ss to b e adjust e d, said co e f 

[00010] Int e gration of (2) r e quir e s an int e gration boundary condition. Advantag e ously, said 

int e gration boundary condition is d e t e rmin e d so that the value D m (h) at th e top of th e cloud cor 
r e sponds to th e pr e d e t e rmin e d valu e for the total numb e r of particl e s at th e top of th e cloud. 

[00011] In an advantag e ous impl e m e ntation, the profil e of th e total numb e r of particl e s n ^fe) 

is d e t e rmin e d by th e following e quation: 

ft ,(h)-7C.Z(h)/D m (h) 6 
wh e r e x is e qual to 25. 4 10~^t 

[0010] — In anoth e r impl e m e ntation, th e m e t e orological paramet e r N o (h) is d e t e rmin e d by th e 
following e quation: 

wh e r e y is e qual to 102 10~^t 

[0011] In a third impl e m e ntation, th e m e t e orological param e t e r corr e sponding to th e profil e 

of th e ice wat e r content IWC(Ii) (in g/m^] is d e t e rmin e d by th e following e quation: 

IWC(h)-wZ(h)/D ffl (h)* 
wh e r e w is e qual to 1.25 IO^ t 

[0012] In a fourth impl e m e ntation, th e m e t e orological paramet e r corr e sponding to th e profil e 

of th e solid pr e cipitation rat e R(h) (mm/h e quival e nt m e lted) is det e rmined by th e following 
e quation: 

R(h)-r.Z(h)/D ffl (h)^ 
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wh e r e r is equal to 4 .698 10 t 



Brief Description of the Drawings 

[0013] ^h eSelected aspects of the invention will be better understood on reading the follow- 
ing description, given with reference to a non-limiting impl e mentation, example and selected 
graphic representations below: 

■ Figur e 1 shows an e xampl e Fig. 1 is a graph of a vertical profile of Z to be inverted (in 
this example, the isotherm 0°C is at ground level); 

Figur e 2 shows th e Fig. 2 is a graph of a profile of Dm resulting from Z being inverted 
by the aggregation model compared with thea conventional estimator; 

Figur e 3 shows th e Fig. 3 is a graph of profiles of NO and of nT resulting from Z being 
inverted by the aggregation model compared with thea conventional hypothesis and ob- 
servations; 

Figure 4 shows Fig. 4 is a graph of the vertical profile of Z to be inverted; 

Figure 5 shows Fig. 5 is a graph of the profile of Dm resulting from Z being inverted by 
the aggregation model compared with the conventional estimator; 

Figur e 6 shows Fig. 6 is a graph of the profiles of NO and nT resulting from the Z being 
inverted by the aggregation model compared with thea conventional hypothesis and ob- 
servations; 

Figur e 7 shows Fig. 7 is a graph of the profile of R resulting from Z being inverted by 
the aggregation model compared with thea conventional estimator; 

Figur e 8 shows Fig. 8 is a graph of the sensitivity of the retrieval of the equivalent pre- 
cipitation rate (mm/h) to the particle density law p(D) oc D-y; and 



5 



Figur e 9 shows Fig. 9 is a graph of the sensitivity of the retrieval of the vertical Doppler 
velocity to the particle density law p(D) oc D-y. 

Detailed Description 

[0014] This disclosure provides methods of estimating precipitation characteristics and, in a 
particular selected aspect, a precipitation rate for solid precipitation, the method comprising an 
acquisition and processing step comprising acquiring a radar image including at least a vertical 
plane of a precipitation zone, and in processing a vertical profile to deliver digital signals repre- 
sentative of reflectivity in the vertical direction h, an integration step comprising integrating 
signals representative of reflectivity to deliver a signal representative of the profile in the vertical 
plane of a mean particle diameter weighted by the mass of each particle, and a determination step 
comprising determining the concentration of the solid particles on the basis of the signals com- 
puted in the preceding steps. 

[0015] Reference will be made to the drawings below as appropriate. Steps in accordance 
with selected aspects of the method may be seen in one sense by reference to Figs. 1 - 3. 
[0016] In a preferred aspect, the integration step comprises determining the variable ZOi) of 
the radar observable in mm 6 /m 3 as a function of the altitude h on the basis of the radar image, 
and determining the mean diameter Dmfti) of the particles by solving the following equation: 

fD m b-5 (\ 1/Z ^ 

= -0.25ke ff aD 10 -18 Z+ I |D m (2) 

fh ™ l6Zfli J 

where: 

Z is the radar observable to be inverted in mm m" ; 
D m is in meters (m); 
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a and b are coefficients specific to particles of the "aggregate" type; for example, the co- 
efficient a is equal to 35184 and the coefficient b is equal to 3.16: 

/k g- is the coefficient of effectiveness of the aggregation process to be adjusted, said coef- 
ficient kg f r being equal to 0.3, 

Integration of (2) uses an integration boundary condition. Advantageously, the integra- 
tion boundary condition is determined so that the value D m (h) at the top of the cloud corresponds 
to the predetermined value for the total number of particles at the top of the cloud. 

In an advantageous implementation, the profile of the total number of particles nfli) is 
determined by the following equation: 

n I (h) = x.Z(h)/D n i (h) 6 
where x is equal to 25.4 10~ 18 . 

In another implementation, the meteorological parameter N o (h) is determined by the fol- 
lowing equation: 

Np(h) = v.ZfhVD^h) 7 
where y is equal to 102 10" 18 . 

In another aspect, the meteorological parameter corresponding to the profile of the ice 
water content IWCfhK'm g/m 3 ) is determined by the following equation: 

IWCOO = wZQiVDnfli) 3 
where w is equal to 1.25 10~ 12 . 

In yet another aspect, the meteorological parameter corresponding to the profile of the 
solid precipitation rate R(h) (mm/h equivalent melted) is determined by the following equation: 

R(h) = r.Z(h)/Dn,(h) 235 
wherein r is equal to 4.698 10' 10 . 
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[0017] The method applies e ss e ntially to stratiform precipitation. It considers that the-ice- 
forming nuclei are activated only at highly negative temperatures, i.e. a at the top of the cloud. 
The-4eelce crystals formed at high altitude settle and grow as they fall, either by sublimation of 
the ambient saturating water vapor, or by collection and freezing of supermelted cloud water 
droplets, or else-by random aggregation of their collisions with other ice crystals. Of the three 
growth processes, only aggregation changes the ice particle concentration. The characterization 
method is based e ssentially on a simplified description of the aggregation mechanism. Th e in 
v e ntion is bas e d on We use a class of "profiler" methods in-that ^inverts the vertical profile of 
reflectivity measured in the ice , so as to extract therefrom the vertical profile of the solid pre- 
cipitation rate. 

[0018] Th eSelected steps of th e inv e rsion m e thod ar e as folio ws include : 

1 - Particle size distribution expressed in "equivalent melted diameter" is assumed to be 
exponential, i.e.: 

N(D)=N 0 exp(-4D/D m ) (1) 
where N(D) is the concentration of particles per cubic meter (m 3 ) and per diameter range, 
and N 0 and D m are the two parameters that characterize distribution. 

2 - The top h max and the base h min of the layer of solid precipitation are determined; 

a. h max is the maximum altitude of the measured reflectivity profile Z(h). 

b. h m in is either the altitude of the isotherm 0°C if the ground-level temperature is 
positive, or it is ground level if the ground-level temperature is negative. 

3 - The profile of the parameter Dm(h) in the range h max to h min is then determined by 
resolving the following differential equation: 
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fo m b. 5 rn/z ^ 

= -0.25keflaD l(r 18 Z+ I |D m (2) 

fh m \6Zfh J 

where: 

- Z is the radar observable to be inverted in nWW 3 : 

- D m is in meters (m); 

- a and b are coefficients specific to particles of the "aggregate" type, equal re- 
spectively to 35184 and to 3.16 on the basis of the observations of J.D. Locatelli 
and P.V. Hobbs , Fall speeds and masses of solid precipitation particles, J. Geo- 
phvs. Res., 79, 2185 - 2197 (1974 V the subject matter of which is incorporated 
herein by reference ; 

- k e ff is the coefficient of effectiveness of the aggregation process to be adjusted 
(the value k e ff= 0.3 seems correct). 

4 - The integration of (2) takes place from the top, where the boundary condition is ex- 
pressed by fixing the total number of particles n T (or the number of ice-forming nuclei 
activated at the top of the cloud). It is possible to take n T (h max ) = 10 6 m' 3 , which makes it 
possible to express the boundary condition D m (h mw J as: 

Dm(H m J=25.4 10' I8 (Z(h max /nj{h m( J) J/6 (3) 

5 - Once the profile D m (h) from h max to h min has been determined, the profiles of the other 
parameters of interest are computed by the following expressions: 

a. Profile of N 0 : N 0 (h)=102.10" 12 Z(h)/D m (h) 7 

b. Profile of the total number of particles n T (h) {(in in 3 }): N T (h)=102.10" l2 Z(h) 
/D m (h) 6 

c. Profile of the ice water content IWC(h) (in g/m 3 }): 

IWC(h)=l .25 . 1 0' 12 Z(h)/D m (h)3 
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d. Profile of the solid precipitation rate R(h) (mm/h equivalent melted). By using 
the terminal fall velocity determined by Locatelli and Hobbs for aggregates: 
Kvt=107.6 D 0 65 (D in m)}), R(h) is expressed by: 
R(h)=4.698. 1 0~ 10 Z(h)/D m (h) 2 35 

[0019] The following description relates to another veFsien selected aspect of tfeea method of 

processing for determining the precipitation rate. 

[0020] For rainfall, the algorithm used is the ZPHI algorithm which is the subject of Fr e nch 
Patent Documents F R 98/00714 and FR 01/09206 , the subject matter of which is incorporated by 
reference . 

[0021] For solid precipitation, the nov e l algorithm that has just be e n d e v e lop e d is th e subj e ct 
ef-the-fol lowing description applies . 

[0022] This algorithm The method for estimating the precipitation rate for solid precipitation 
b e longs to th e includes a class of "profiler" algorithms ifl-that it-inverts the vertical profile of 
reflectivity measured in the ice , so as to extract therefrom the vertical profile of the solid precip- 
itation rate. 

[0023] The algorithm method applies e ss e ntially to stratiform precipitation. It considers that 
the ice-forming nuclei are activated only at highly negative temperatures, i.e. a at the top of the 
cloud. The ice crystals formed at high altitude settle and grow as they fall, either by sublimation 
of the ambient saturating water vapor, or by collection and freezing of supermelted cloud water 
droplets, or else by random aggregation of their collisions with other ice crystals. Of the three 
growth processes, only aggregation changes the ice particle concentration. The algorithm is 
based essentially on a simplified description of the aggregation mechanism. : Fh eS elected steps 
of the inversion method are as follows: 
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Particle size distribution expressed in "equivalent melted diameter" is assumed to 
be exponential, i.e.: 

NfD^NtfixpHD/DnJ (1) 
where D is the equivalent melted diameter of the ice particle; 
N(D) is the concentration of particles per m 3 and per diameter range; and 
N 0 and D m are the two parameters that characterize distribution. 

The top hmax and the base h m i n of the layer of solid precipitation are determined; 
h max is the maximum altitude of the measured reflectivity profile Z(h). 
hmin is either the altitude of the isotherm 0°C if the ground-level temperature is 
positive, or it is ground level if the ground-level temperature is negative. 
[0024] The profile of the parameter D m (h) in the range h max to h min is then determined by re- 
solving the following differential equation, whose solution can be determined analytically: 

dOm 4 d ke ff J(b,d)a W) 1 dZ 

= -24 ZD + D m (2) 

ah T(4+d) a(a-l) m (a-l)Zdh 

where: 

foo fx r i 

J(b,d)= I exp(-4x)x b I I x d -x' d lexpHx^dx'dx 

Jo Jo L J 

describes the collision frequency; 

Z is the radar observable to be inverted in WW 3 : 
D m is in m; 

- a, b, c, and d are coefficients dependent on the density law F(D) for ice particles, 
which law is assumed to vary at Z)" 7 (where 7 can take values in the range 0.25 to 1.1 
depending on the type of particles). These coefficients are drawn from Mitchell's 
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Theory (Journal of Atmospheric Sciences, 53, 12, 1996) . the subject matter of which 
is incorporated by reference, for representing the following power laws: 

v t (D) = cD d (where v, is the terminal fall velocity of the particle of equivalent 
melted diameter D): 

A = aD b (where A is the effective cross-section of the particle of equiv- 
alent melted diameter D); 
With coefficients depending on frequency and enas a function of (D), represent- 
ing by a power law the relationship between reflectivity (measured parameter), the con- 
centration of the particles, and their mean equivalent diameter: 
Z = aN 0 D m d 

kefr is the coefficient of effectiveness of the aggregation process to be adjusted 
(the value 1^= 0.3 seems correct immediately above the isotherm 0°C). 
[0025] The integration of (2) takes place from the top, where the boundary condition is ex- 
pressed by fixing the total number of particles n T (or the number of ice-forming nuclei activated 
at the top of the cloud). It is possible to take nj (h max ) = 10 6 m" 3 , which makes it possible to ex- 
press the boundary condition D m (h max ) as: 

1 

ri Z(h max ) 1 a-\ 
D m (h max )= I | (3) 

[0026] Once the profile D m (h) from h max to h m j n has been determined, the profiles of the other 
parameters of interest are computed by the following expressions: 
Profile of N 0 : 
Z(h) 

N 0 (h)= 

dD m (h) & 
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Profile of the total number of particles n T (h) {{in m' 3 }): m{h)=025No(h)D m (h) 
Profile of the ice water content IWC(h) (in g/m 3 }}: 

IWC(h)={\2212)\tfN 0 (h)D m (hf 
Profile of the solid precipitation rate R(h) (mm/h equivalent melted): 

T(4+d) 

R(h)=(lM5)l0 6 c N 0 (h)D m (h) (4+d) 

4 (4+d) 

Vertical profile of the Doppler velocity V D (h) (in meters per second (m/s)), given 
by a power law dependent on density _(D): 
V D (h)=pD m (hf 

[0027] Because of the dependency of the coefficients relative to the density law F(D) 9 the 
retrieval of the precipitation rate depends critically on the parameter y, as shown by figw eFig. 7. 
The density law thus constitutes a key parameter. 

[0028] By operating in vertical firing, it is possible to measure the Doppler velocity V D (h), 
and te-compare it with V D (h) computed by the algorithm. The Doppler velocity profile is dis- 
criminating, as shown by FigweFig. 8. By successive adjustments, it is thus possible to deter- 
mine the parameter % of the law F(D) for which consistency between V D (h) as computed and 
Vo(h) as measured is achieved. Figur eFig. 4 shows an example of a profile inverted by the 
algorithm (in thisthat example, the isotherm 0°C is at ground level). Figuf eFig. 5 shows the 
curves resulting from application of the aggregation model and of th ea conventional model, 
making it possible for the profile of D M resulting from the inversion of Z by the aggregation 
model to be compared with the conventional estimator. 

[0029] Figur eFig. 6 shows the-curves of the profiles as a function of altitude (m) above the 
isotherm 0°C. It enables the profiles of N 0 and n T resulting from Z being inverted by the aggre- 
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gation model to be compared with the conventional hypothesis and observations. FigweFig. 7 
shows the profile of R resulting from Z being inverted by the aggregation model compared with 
the conventional estimator. FigweFig. 8 shows the sensitivity of the retrieval of the equivalent 
participation rate (mm/h) to the particle density law p(D) oc D" Y . FigweFig. 9 shows the sen- 
sitivity of the retrieval of the vertical Doppler velocity to the particle density law p(D) oc D" Y . 
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